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We review discoveries in the nonlinear dynamics of curved spacetime, largely made possible by 
numerical solutions of Einstein's equations. We discuss critical phenomena and self-similarity in 
gravitational collapse, the behavior of spacetime curvature near singularities, the instability of black 
strings in 5 spacetime dimensions, and the collision of four-dimensional black holes. We also discuss 
the prospects for further discoveries in geometrodynamics via observation of gravitational waves. 


I. INTRODUCTION 


In the 1950s and 60s, John Archibald Wheeler [1] spec- 
ulated that curved, empty spacetime could exhibit rich, 
nonlinear dynamics — which he called geometrodynam- 
ics — analogous to the writhing surface of the ocean in a 
storm. Wheeler exhorted his students and colleagues to 
explore geometrodynamics by solving Einstein's general 
relativistic field equations. 

In 1965, Yakov Borisovich Zel'dovich, with his young 
protégés Andrei Doroshkevich and Igor Novikov [2], gave 
strong evidence that, when a highly deformed star col- 
lapses to form what would later be called a black hole, 
the hole's curved spacetime, by its nonlinear dynamics, 
will somehow shake off all the deformations, thereby be- 
coming a completely smooth, axially symmetric hole. 

The Wheeler/Zel'dovich challenge of exploring ge- 
ometrodynamics in general, and for black holes in par- 
ticular, has largely resisted analytic techniques. In the 
1980s and 90s, that resistance motivated the authors and 
our colleagues to formulate a two-pronged attack on ge- 
ometrodynamics: numerical solutions of Einstein’s equa- 
tions to discover general relativity's predictions, and ob- 
servations of gravitational waves from black-hole births 
and black-hole collisions to test those predictions. The 
numerical simulations have now reached fruition, and the 
observations will do so in the near future. 

In this article—dedicated to the memory of Zel'dovich 
and Wheeler (who deeply respected each other despite 
cold-war barriers, and who were primary mentors for one 
of us, Kip Thorne)—we shall present an overview of some 
of the most interesting things we have learned about ge- 
ometrodynamics from numerical simulations, and a pre- 
view of what gravitational-wave observations may bring. 
More specifically: 

We shall describe geometrodynamic discoveries in four 
venues: gravitational implosion, where a phase transi- 
tion, discrete self-similarity and critical behavior have 
been observed (Sec. mp; the dynamics of spacetime near 
singularities, where richly chaotic behavior has been ob- 
served (Sec. mm); the unstable evolution of a black string 
in five spacetime dimensions, where a dynamical se- 
quence of strings linking black holes has been observed 
(Sec. [IV]; and collisions of black holes in four spacetime 
dimensions, where dynamical interactions of tidal ten- 
dexes and frame-drag vortexes have been observed (Sec. 


[V). Then we shall briefly describe the prospects for ob- 
serving some of these phenomena via gravitational waves 


(Sec. [V1]. 


II. GRAVITATIONAL COLLAPSE: PHASE 
TRANSITION AND CRITICAL BEHAVIOR 


The first numerical simulations to exhibit rich ge- 
ometrodynamics were in 1993, by Matthew Choptuik 
[3], who was then a postdoc at the University of Texas. 
Choptuik simulated the spherical implosion (Fig. of 
a linear, massless scalar field (satisfying OY = 0). The 
field’s energy, momentum and stress (which are quadratic 
in the field) generated spacetime curvature, with which 
the field interacted via the curvature’s influence on its 
wave operator []. That interaction produced surprising 
nonlinear dynamics: 


FIG. 1: The implosion of a scalar wave V with amplitude p 
and a particular waveform. 


If the wave amplitude p, for some chosen ingoing wave- 
form, was larger than some critical value p,, the implo- 
sion produced a black hole. If p was smaller than p., the 
imploding waves passed through themselves, travelled 
back outward, and dispersed. For p — p., the imploding 
waves interacted with themselves nonlinearly (via their 
spacetime curvature), producing a sequence of frequency 
doublings and wavelength halvings, with an intriguing, 
discretely self-similar pattern that was independent of 
the initial, ingoing waveform. Waves with ever decreasing 
wavelength emerged from the “nonlinearly boiling" field, 
carrying away its energy, and ultimately leaving behind 
what appeared to be an infinitesimal naked singularity (a 
region with infinite spacetime curvature, not hidden by 


8 black-hole horizon). One year after Choptuik's simu- 
lations, the mathematician Demetrios Christodoulou 
developed a rigorous proof that the endpoint, for p = px, 
was, indeed, a naked singularity. 

The transition of the implosion's endpoint, from black 
hole for p > p,, to naked singularity for p = p,, to wave 
dispersal for p « p., was a phase transition analogous 
to those that occur in condensed-matter physics. And, 
as in condensed matter, so also here, the phase transi- 
tion exhibited scaling: For p slightly larger than p,, the 
mass of the final black hole scaled as Mpu x (p — p.)^. 
For p slightly below p., the radius of curvature of space- 
time at the center of the ^boiling region" reached a 
minimum value, before field dispersal, that scaled as 
Rmin = (R Ruvop)mak’ x (p. — p)?, with the same 
numerically measured exponent 8 = 0.374. Here Ryuvep 
is the Riemann curvature tensor. 

Choptuik's discovery triggered many follow-on simula- 
tions. Most interesting to us was one by Andrew Abra- 
hams and Charles Evans [5], later repeated with higher 
resolution by Evgeny Sorkin [6]. In their simulations, 
the imploding, spherically symmetric scalar field was re- 
placed by an imploding, axially symmetric (quadrupolar) 
gravitational wave — so they were dealing with pure, vac- 
uum spacetime as envisioned by Wheeler. Once again, 
there was a critical wave amplitude p.; and for p near 
p. the behavior was similar to the scalar-wave case, to 
within numerical error: for p = p,, strong evidence for 
discretely self-similar evolution leading to an infinitesi- 
mal final singularity; for p > p., the same black-hole 
mass scaling Mpu « (p — p.)^; for p < px, the same 
spacetime-curvature scaling R œ (p, — p)?; and to within 
numerical accuracy, the scaling exponent was the same: 
6 = 0.38 for the quadrupolar gravitational waves, and 
B = 0.374 for the spherical scalar wave. This is remi- 
niscent, of course, of the universality one encounters in 
condensed-matter phase transitions. 

For a detailed review of these and many other studies of 
critical gravitational implosion, see an article by Carsten 
Gundlach [7]. 


III. GENERIC SPACE-TIME SINGULARITIES 
A. BKL Singularity 


The singular endpoint of the implosions described 
above is non-generic, in the sense that no singularity oc- 
curs if p, is only infinitesimally different from p. 

However, there are other situations that lead to generic 
singularities.! This was demonstrated analytically in the 


1 Perhaps the most important generic singularity for astrophysics 
is one that arises when matter with negligible pressure is present. 
Leonid Petrovich Grishchuk in 1967 [8| showed that the matter, 
generically, undergoes gravitational collapse to form flat pan- 


1960s by Roger Penrose, Stephen Hawking and others, 
using tools from differential topology [10]. In 1969-70, 
Vladimir Belinsky, Isaac Khalatnikov and Evgeny Lif- 
shitz [II] used approximate differential-geometry tech- 
niques to deduce the geometrodynamical behavior of 
spacetime as it nears one generic type of singularity, a 
type now called BKL. 

In the 1970s, 80s and 90s, there was much skepticism 
in the US and UK about this BKL analysis, because its 
rigor was much lower than that of the Penrose-Hawking 
singularity theorems. (This lower rigor was inevitable, 
because the geometrodynamical approach to a singular- 
ity is very complex—see below—and deducing its details 
is much more difficult than proving that a singularity oc- 
curs.) As a result, the BKL geometrodynamics came to 
be called, in the West, the BKL conjecture. 

There was little hope for proving or disproving this 
"conjecture" by analytical techniques, so the skeptics 
turned to numerical simulations for probing it. After 
nearly a decade of code development, David Garfinkle 
in 2003 [12] carried out simulations which demonstrated 
that Belinsky, Khalatnikov and Lifshitz had been cor- 
rect. The geometrodynamical evolution, approaching a 
BKL singularity, is as they predicted, though with one 
additional feature that they had missed: a set of nonlocal 
spikes in the spacetime curvature E 

The BKL geometrodynamics can be described in terms 
of tidal-gravity observations by observers who fall into 
the BKL singularity on timelike geodesics; Fig.[2] As two 
observers, A and B, approach the singularity, they lose 
causal contact, in the sense that, after passing through 
A's particle horizon (at point P in the diagram), B can 
no longer influence A. This causal decoupling is so strong, 
in the BKL spacetime, that spatial derivatives cease hav- 


Singularity 


FIG. 2: The world lines of two observers, A and B, falling into 
a BKL singularity (the solid lines), and the particle horizon 
of observer A (dashed lines; a past light cone). Events outside 
the particle horizon can never influence observer A. 


cakes with infinite density and curvature; and in 1970 Zel'dovich 
[9] showed that, astrophysically, pressure halts the collapse before 
the infinities but after the pancake structure has been strongly 
established. A few years later Zel'dovich realized that these pan- 
cakes, seen edge on, are filamentary structures that astronomers 
observe in the distribution of galaxies on the sky. 


ing significant influence on the geometrodynamics, as the 
singularity is approached—there is a spatial decoupling— 
and as a result, it turns out, there is no correlation be- 
tween the late-time observations of adjacent observers. 


Each observer’s experience, when approaching the sin- 
gularity, can be described in terms of the tidal gravita- 
tional field £j; that he feels. This tidal gravitational field 
has components, in the observer’s local Lorentz frame, 
that are equal to the space-time-space-time part of the 
Riemann curvature tensor. 


Ejk = Ryoko - (1) 


The physical manifestation of the tidal field is a relative 
gravitational acceleration 


Aa; = — Gk AL (2) 


of particles with vector separation Azp. (The tidal field 
acquires its name from the fact that the Sun’s and Moon’s 
tidal field produces the tides on the Earth’s oceans. In 
the Newtonian limit, the tidal field is £j = 07/02; 02x, 
where ® is the Newtonian gravitational potential.) 


Being a symmetric tensor, the tidal field can be de- 
scribed by three orthogonal principal axes (unit vectors 
ei, e5, €3), and their eigenvalues, £;; = e; - E. e;. If 
€i; < 0, then an object is tidally stretched along its 
e; principal axis, and similarly for the other principal 
axes. If £;; > 0, the object is tidally squeezed along ez. 
'The tidal field, in vacuum, is trace-free, so its eigenvalues 
must sum to zero, which means there must be a squeeze 
along at least one principal axis and a stretch along at 
least one. 


Figure |3| shows the pattern of stretches and squeezes 
experienced by an observer when falling into the BKL 
singularity. The pattern is divided, in time, into cycles 
and eras. During a single cycle, there is a stretch along 
one axis and a squeeze along the other two. Between 
cycles, the stretch axis switches to squeeze and the more 
strongly squeezing axis switches to stretch. At the end of 
each era, the axes rotate in some manner, relative to the 
observer's local Lorentz frame, and the pattern begins 
over. The number of cycles in each era and the details of 
their dynamics are governed by a continued-fraction map 
that is chaotic, in the technical sense of chaos (extreme 
sensitivity to initial conditions). This chaos plays a key 
role in destroying correlations between the observations 
of adjacent observers as they approach the singularity. 


The full details of this, as worked out by Belinsky, 
Khalatnikov and Lifshitz [II], occasionally are violated 
(numerical simulations reveal): A cycle gets skipped, and 
during that skip, there is an extreme spike in the tidal 
field that is more sensitive to spatial derivatives than 
expected, and whose details are not yet fully understood. 
See [13] and references therein. 


STRETCH 
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FIG. 3: The qualitative pattern of tidal stretches and squeezes 
experienced by an observer when falling into a BKL singular- 
ity. The three eigenvalues of the tidal field are plotted verti- 
cally as functions of time, with one axis shown solid, another 
dashed, and the third dotted. Adapted from 


B. Singularities Inside a Black Hole 


This BKL behavior is speculated to occur in the core 
of a young black hole. However, numerical simulations 
are needed to confirm or refute this speculation. 

As the black hole ages, the singularity in its core is 
speculated to break up into three singularities, one of 
BKL type, and two that are much more more gentle 
than BKL. This speculation arises from the expecta- 
tion that, just as the hole’s exterior spacetime geometry 
settles down into the quiescent, axially symmetric state 
described by the Kerr metric, so its interior will settle 
down into the Kerr-metric state, except near two special 
regions called Cauchy horizons, where the Kerr metric 
is dynamically unstable. Nonlinear geometrodynamics 
is thought to convert those Cauchy horizons into null, 
generic singularities (Fig. (4). 

These singularities are null in the sense that ingoing or 
outgoing photons, in principle, could skim along them, 
not getting captured. The ingoing singularity (called a 
mass inflation singularity) is generated, according to ap- 
proximate analytical analyses [I6] [I7], by radiation and 
material that fall into the black hole and pile up along 
the ingoing Cauchy horizon, there gravitating intensely. 
The outgoing singularity (called a shock singularity [15]) 
is generated by ingoing radiation that backscatters off 
the hole’s spacetime curvature, and then travels outward, 
piling up along the outgoing Cauchy horizon, there grav- 
itating intensely. In both cases, the piling-up stuff could 


FIG. 4: Penrose diagram, depicting the causal structure out- 
side and inside an old black hole, as best we understand it 
today. Ingoing and outgoing null rays (hypothetical photons) 
travel along 45 degree lines, and a conformal transformation 
has been used to compress our universe into a finite diamond 
in the diagram. The Kerr spacetime is shaded. The true 
spacetime is superposed on the Kerr spacetime; it is the re- 
gion bounded by the center of the imploding star (thin left 
line), the BKL singularity (thick horizontal line), the mass- 
inflation singularity (dashed line), and the infinities of our 
universe: future timelike infinity labeled J, future null infin- 
ity labeled Z,, spacelike infinity labeled Io, past null infinity 
labeled Z., and past timelike infinity labeled J_. It might 
well be that the true spacetime ends at the shock singularity, 
and there is no BKL singularity beyond [15]. 


be gravitational waves rather than material or nongravi- 
tational radiation, in which case we are dealing with pure 
vacuum spacetime curvature—pure geometrodynamics. 

Neither of these null singularities is oscillatory, but at 
both of them, the spacetime curvature goes to infinity 
(the radius of curvature of spacetime R goes to zero). 
'This divergence of curvature happens so quickly at the 
shock singularity, that objects might be able to pass 
through it, though with a destructive net compression 
along two dimensions and net stretch along the third. If 
so, they will likely then be destroyed in the BKL singu- 
larity. 

'The mass-inflation singularity is expected also to pro- 
duce only a finite, not infinite, net compression and 
stretch of objects falling through. If anything survives, 
its subsequent fate is totally unknown. 

These (highly informed) speculations, which arise from 
extensive, approximate analytically analyses, mostly per- 
turbation theory, will be tested, in the coming few years, 
by numerical simulations. And just as the BKL con- 
jecture missed an important phenomenon (the curvature 


spikes), so these speculations about the geometrodynam- 
ics of black-hole interiors may fail in some modest, or 
even spectacular way. 

For far greater detail on what we now know and spec- 
ulate about the interiors of black holes and the bases for 
that knowledge and speculation, see and references 
therein. 


IV. BLACK STRING IN FIVE SPACETIME 
DIMENSIONS 


A remarkable example of geometrodynamics has been 
discovered by Luis Lehner and Frans Pretorius (18|, in 
numerical simulations carried out in five spacetime di- 
mensions. 

Lehner and Pretorius began with a black string in its 
equilibrium state. This black string is a vacuum solution 
of Einstein's equations in five spacetime dimensions, with 


metric 
2M dr? 
ds? = —(1 t? 
° ( r E 1—2M/r 
+r?(d6? + sin? 0d?) + dz? . (3) 
This is precisely a four-spacetime-dimensional 


Schwarzschild black hole, translated along the z 
axis in the fifth (spatial) dimension. The event horizon 
is at r = 2M; at fixed time t, it is a cylinder with 
spherical cross section, R x S?. 

In 1993, Ruth Gregory and Raymond LaFlamme [19] 
proved, analytically, that such a black string is unstable 
against linear, axisymmetric perturbations with wave- 
length longer than about 1.2 times the string's circum- 
ference. But little definitive was known about the in- 
stability's nonlinear, geometrodynamical evolution until 
Lehner and Pretorius's 2010 simulations |18]. They re- 
vealed that the string's horizon evolves as depicted in 
Fig. 

The string develops a sausage instability—analogous 
to that of a magnetically confined plasma in a Z-pinch, 
and the Rayleigh-Plateau instability for a cylinder of 
fluid confined by its surface tension—but with outgo- 
ing gravitational waves carrying off energy. This insta- 
bility gives rise to a chain of five-spacetime-dimensional 
black holes linked by segments of shrunken black string— 
segments whose circumferences are far smaller than that 
of the original string. The instability then repeats 
on each shrunken string, producing smaller black holes 
linked by segments of more extremely shrunken string. 
Each successive sausage instability produces its smaller 
black holes on an evolution timescale proportional to 
the string's circumference. These successively shorter 
timescales converge: An infinite sequence of instabilities, 
in finite time, presumably ends in a naked singularity. 

Because these simulations assumed 2-sphere symmetry 
from the outset, we cannot be certain that their predic- 
tions are the black string's true geometrodynamical evo- 
lution. To learn the true evolution, we need simulations 
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FIG. 5: (Color online.) A sequence of snapshots from a simulation of the geometrodynamical evolution of a black string in 


5 spacetime dimensions, by Lehner and Pretorius |18]. 


Each snapshot is an embedding diagram of the black string’s event 


horizon: the horizon’s intrinsic geometry is the same as the intrinsic geometry of the depicted surface in flat space. 


in five spacetime dimensions, that do not assume any 
symmetry. Such simulations are beyond current capabil- 
ities, but may be possible in a decade or so. In the mean- 
time, it is conjectured that the Lehner-Pretorius evolu- 
tion (Fig. 5) is the true evolution, because black strings 
have been proved stable against all linear nonspherical 
perturbations. 


V. BLACK-HOLE COLLISIONS 


Recent advances in numerical simulations have enabled 
the study of geometrodynamics in the violent collisions 
of two black holes—including the generation of gravi- 
tational waves, and the relaxation of the remnant to a 
single, quiescent, spinning black hole (as predicted by 
Zel'dovich, Doroshkevich and Novikov [2]). We and col- 
laborators have developed a new set of so-called vor- 
tez/Lendex tools for visualizing this black-hole geometro- 
dynamics [20]. We will first describe these tools, and 
then we will use them to visualize the geometrodynamics 
of black-hole collisions. 


A.  Vortex-tendex tools 


The gravitational field felt by a local observer is de- 
scribed by the Riemann curvature tensor Ryvop. Any 
such observer, freely falling or accelerated, can decom- 
pose Riemann into an “electric” part, Ejk, defined by 
Eq. {i}, and a “magnetic” part, Bjk, defined by 


Bik = €jpqRpqko - (4) 


Here the indices are components on the observer’s local, 
orthonormal basis, the index 0 refers to the time com- 
ponent (i.e., the component along the observer’s world 
line), Latin indices refer to the observer’s three spatial 
components, and €jpq is the spatial Levi-Civita tensor. 


Both £j; and Bj, are symmetric and trace free in vac- 
uum (the situation of interest to us). 

As discussed in Section £j is called the tidal 
field, and describes the tidal stretching and squeezing 
experienced by objects in the observer's local reference 
frame, according to Eq. (2). The “magnetic” quantity 
Bj, is called the frame-drag field. The physical manifes- 
tation of this field is a relative precession, or dragging 
of inertial frames: two gyroscopes with vector separa- 
tion Az, will precess relative to each other with angular 
velocity 


AQ; = BjpAxp. (5) 


This differential frame dragging is a general relativistic 
effect with no analogue in Newtonian gravity. Its global 
(non-differential) analog is one of the two relativistic ef- 
fects recently measured by Gravity Probe B [21]. 

Note that the decomposition of the Riemann tensor 
into Ejk and Bj, depends on the observer’s reference 
frame. Different observers at the same location, moving 
relative to each other, will measure different tidal fields 
and frame-drag fields. This is the same situation as for 
classical electromagnetism, in which the electromagnetic 
tensor F y can be decomposed into the familiar electric 
and magnetic field vectors, E; = Fjo and Bj = €jpqF pq; 
in the same observer-dependent manner. This correspon- 
dence between gravitation and electromagnetism moti- 
vates the use of the names “electric” and “magnetic” for 
Ejk and Bjk. 

The frame-drag field Bj; and the tidal field £j; are 
useful in describing geometrodynamics at the surface of 
a black hole (its event horizon). If n is a unit vector 
normal to the hole's horizon, with spatial components 
n', then we define the horizon tendicity Enn = nin*E;x, 
as the normal-normal component of the tidal field. For 
positive horizon tendicity, an object is tidally squeezed 
normal to the horizon, and for negative horizon tendicity, 
an object is tidally stretched normal to the horizon. This 
is illustrated in the left panel of Fig. [e] We call a region 
on the horizon with large tendicity a horizon tendez. 
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FIG. 6: (Color online.) A spinning black hole. Left: Colors 


represent horizon tendicity Enn. There is a positive (blue 
or light gray) tendex on each of the poles, and a negative 
(red or dark gray) tendex on the equator. Other (green or 
medium gray) regions have small tendicity. Right: Colors 
represent horizon vorticity Bnn. There is a negative (red or 
dark gray) vortex on the north pole, and a positive (blue or 
light gray) vortex on the south pole. Other (green or medium 
gray) regions have small vorticity. The spin vector points out 
of the north pole. 


We can similarly define the horizon vorticity Bnn = 
nj në Bjk as the normal-normal component of the frame- 
drag field at the surface of a black hole. For positive hori- 
zon vorticity, an object experiences a clockwise twist; for 
negative horizon vorticity, the twist is counterclockwise. 
We call a region on the horizon with large vorticity a 
horizon vortex. The horizon vorticity of a spinning black 
hole is illustrated in the right panel of Fig. [6] 

We now turn to vortex/tendex tools in regions away 
from the horizon. In Section we discussed how, 
at any point, Ejk can be described by three orthogonal 
eigenvectors (unit vectors ej, e5, e4), and their eigenval- 
ues, £5; = e5 : E: ej. We call the eigenvalue €;; the ten- 
dicity associated with the corresponding eigenvector e;; 
the tendicity measures the tidal stretching or squeezing of 
an object oriented along the eigenvector. In analogy with 
electric field lines, we define tendex lines as the integral 
curves along each of the three eigenvectors e;. Whereas 
in electromagnetism there is a single electric field line 
passing through each point, in geometrodynamics there 
are generically three tendex lines passing through each 
point: one tendex line for each of the three eigenvectors 
of Ejk. Because (in vacuum) £j; has vanishing trace, the 
eigenvalues must sum to zero, so generically there exist 
both positive and negative tendex lines at each point. 

The left panel of Fig. [7| shows the tendex lines outside 
a rapidly rotating black hole. A collection of tendex lines 
with particularly large tendicity is called a tendez. The 
rotating hole has a fan-shaped, stretching (red or dark 
gray) tendex sticking out of its equatorial horizon tendex, 
and a poloidal squeezing (blue or light gray) tendex that 
emerges from its north polar horizon tendex, and reaches 
around the hole to its south polar horizon tendex. 

As with Ejk, the frame-drag field Bj, can be described 
by three orthogonal eigenvectors and their eigenvalues. 
An integral curve of one of these eigenvectors is called a 
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FIG. 7: (Color online.) Tendex lines (left panel) and vortex 
lines (right panel) near a spinning black hole. Lines with pos- 
itive eigenvalues (tidal squeeze or clockwise twist) are shown 
blue or light gray, and lines with negative eigenvalues (tidal 
squeeze or counterclockwise twist) are shown red or dark gray. 
At each point in space there are three intersecting tendex lines 
and three intersecting vortex lines. 


vortex line, and the corresponding eigenvalue is the vor- 
ticity associated with that vortex line. The vorticity of a 
vortex line describes the twist, or differential frame drag- 
ging, experienced by an object oriented along that vortex 
line: positive vorticity corresponds to a clockwise twist, 
and negative vorticity corresponds to a counterclockwise 
twist. 

The vortex lines associated with a spinning black hole 
are shown in the right panel of Fig. [7] A counterclockwise 
(red or dark gray) vortex (collection of large-negative- 
vorticity lines) emerges from the horizon’s north polar 
vortex, reaches around the south polar region and de- 
scends back into the horizon’s north polar vortex. Simi- 
larly, a clockwise (blue or light gray) vortex emerges from 
the south polar horizon vortex, reaches around the north 
polar region and descends back into the south horizon 
vortex. 

The black holes of Figs. [6| and [7| have stationary 
(time-independent) vortex and tendex structures. Vor- 
tex and tendex lines, and their associated vortexes and 
tendexes, can also behave dynamically. The equations of 
motion for £;; and Dj, are similar to Maxwell’s equa- 
tions. Like Maxwell, they have wavelike solutions— 
gravitational waves—in which energy is fed back and 
forth between Ejk and Bj,. Figure g shows vortex and 
tendex lines for a plane gravitational wave without any 
nearby sources. As the wave passes an observer, the ten- 
dicities and vorticities oscillate in sign with one oscilla- 
tion per gravitational wave period, leading to an oscilla- 
tory stretch and squeeze in horizontal and vertical direc- 
tions, and an oscillatory twist in diagonal directions, out 
of phase with the stretch and squeeze. 


B.  Black-hole collisions 


We shall illustrate the geometrodynamical richness of 
black-hole collisions by describing the vortex/tendex be- 


FIG. 8: (Color online.) Snapshot of vortex and tendex lines 
for a gravitational plane wave traveling into the page. Two 
orthogonal sets of tendex lines (left) are oriented 45 degrees 
with respect to two orthogonal sets of vortex lines (right). 
The third set of vortex (and tendex) lines is normal to the 
page, with zero vorticity (tendicity). 


haviors in several numerical simulations. All these simu- 
lations were performed by members of the Collaboration 
to Simulate Extreme Spacetimes (SXS), which included 
numerical relativists from Caltech, Cornell, the Canadian 
Institute for Theoretical Astrophysics, and Washington 
State University at the time of these simulations, and has 
since been expanded. We are members of this collabora- 
tion, and one of us (Scheel) played a significant role in 
most of these simulations. 


1. Head-on collision of two black holes with transverse spins 


Our first example is a simulation [22] of the head-on 
collision of two transversely spinning black holes, de- 
picted in Fig. [9] 

As the black holes merge, the vortexes retain their in- 
dividuality. When the four vortexes (one pair from each 
hole) discover each other, all attached to the same hori- 
zon, they begin to interact: each one tries to convert 
those adjacent to it into a replica of itself. As a result, 
they exchange vorticities; each oscillates back and forth 
between clockwise and counterclockwise. At the moment 
when all are switching direction, so the horizon momen- 
tarily has zero vorticity, the vortex lines have popped 
off the horizon and joined onto each other, creating a 
toroidal structure, much like a smoke ring, that is begin- 
ning to travel outward. Simultaneously, adjacent to the 
horizon much of the oscillation energy is stored in ten- 
dexes, which then regenerate the horizon vortexes, but 
with twist directions reversed. As the toroidal bundle of 
vortex lines travels outward, its motion generates toroidal 
tendex lines, intertwined with the vortex lines in just such 
a manner as to become, locally, the gravitational-wave 
structure described in Fig. [8] above. 

The process repeats over and over, with successive, 
toroidal, tendex/vortex structures being ejected and 
morphing into gravitational waves. The waves carry away 


oscillation energy, and some oscillation energy goes down 
the hole, so the oscillations die out, with an exponentia- 
tion time of order an oscillation period. 


2. Collision of identical, spinning black holes, in an 
inspiraling circular orbit 


For collisions of orbiting (i.e. non head-on) black holes, 
the vortex and tendex lines similarly travel to the wave 
zone and become gravitational waves. 

The left panel of Fig. shows a schematic diagram 
of the horizon vortexes and the vortex lines for the colli- 
sion of two orbiting, spinning black holes that are about 
to merge. Just after merger, the horizon vortexes retain 
their individuality, and travel around the horizon of the 
merged black hole, trailing their vortex lines outward and 
backward in a pattern like water from a spinning sprin- 
kler head (shown schematically on the right of Fig. [10]. 
In the wave zone, the vortex lines acquire tendex lines 
and become gravitational waves. 

Similarly, the near-field tendex lines, attached to the 
merged hole's horizon tendexes, trail backward and out- 
ward in a spiral pattern, acquiring accompanying vortex 
lines as they travel, and becoming gravitational waves. 
These horizon-tendex-generated gravitational waves have 
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FIG. 9: (Color online). The event horizons and vortex lines of 
spinning black holes, colliding head-on with transverse spins; 
from the simulation reported in [22]. (a) Horizon vortexes and 
spin directions just before merger. (b) Horizon vortexes just 
after merger, which retain their individuality. (c) Vortex lines 
linking horizon vortexes of same polarity (red to red or dark 
gray to dark gray; blue to blue or light gray to light gray). 
Lines are color coded by vorticity (different scale from hori- 
zon). (d) Sloshing of near-zone vortexes generates toroidal 
vortex loops (two shown here) composed of orthogonal vor- 
tex lines, traveling outward as gravitational waves; these are 
accompanied by intertwined tendex lines (not shown). The 
horizon, with attached vortex lines, is visible in the center. 
Figure adapted from [20]. See for a movie of horizon 
vortexes for this simulation. 


the opposite parity from the horizon-vortex-generated 
waves; and there is a remarkable duality between the 
two sets of waves [24] 

Figure [I0] is schematic. For a more precise depiction, 
we focus on the merged hole at late times, when it is 
weakly perturbed from its final, Kerr-metric state, and 
its perturbations are predominantly / — 2, m — 2 quasi- 
normal modes [24]. Then the perturbations of the frame- 
drag field, 0B;;, generated by the horizon vortexes, have 
the vortex lines and vorticities shown in the left panel of 
Fig. and the perturbations of the tendex field, 0£;;, 
generated by the horizon tendexes, have the tendex lines 
and tendicities show in the right panel. Notice that the 
two figures are very nearly identical, aside from sign (in- 
terchange of red and blue, or dark and light gray). This 
is due to the (near) duality between the two sets of per- 
turbations. 


3.  Extreme-kick black-hole collision 


An interesting example of geometrodynamics, and of 
the interplay between vortexes and tendexes, is the 
"extreme-kick" black-hole collision first simulated, not 
by our SXS collaboration, but by Campanelli et al. [25] 
and others [26] [27]. Our collaboration has repeated their 
simulations, in order to extract the vortex and tendex 
structures. 

In these simulations, two identical black holes merge 
from an initially circular orbit, with oppositely directed 
spins lying in the orbital (x,y) plane. The name 
“extreme-kick” arises because gravitational waves gen- 
erated during the merger carry linear momentum prefer- 
entially in either the +z or —z direction, resulting in a 
gravitational recoil of the remnant black hole with speeds 
as high as thousands of km/s. 'The magnitude and direc- 
tion of the recoil depends on the angle between holes' 
identical spin axes, and the holes separations, at the mo- 
ment of merger. This angle can be fine-tuned (for in- 
stance, to produce maximum recoil in the +z direction) 
by adjusting the initial conditions in the simulation. 

To understand the mechanism of the recoil, consider 


FIG. 10: (Color online.) Left: Schematic of vortex lines 
and horizon vortexes for two orbiting, spinning black holes 
about to merge. Right: Schematic of vortex lines of the 
remnant black hole just after merger, showing vortex lines 
extending to large distances; the entire pattern is rotating 
counterclockwise. 
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FIG. 11: (Color online.) Vortex and tendex structures for 
deviations from the final, Kerr black hole, at late times after 
the merger of a black-hole binary. These are the structures in 
the final hole's equatorial plane, and the final hole has a di- 
mensionless spin S/M? = a/M = 0.945. Left: Perturbations 
generated by horizon vortexes—vortex lines (solid black for 
clockwise, dashed for counterclockwise), and vorticity of the 
dominant vortex at each point (colored blue or light gray for 
clockwise and red or dark gray for counterclockwise). Right: 
Perturbations generated by horizon tendexes—tendex lines 
(solid black for squeeze, dashed for stretch), and tendicity of 
the dominant tendex (blue or light gray for squeeze and red 
or dark gray for stretch). Adapted from Fig. 12 of [24]. 


the remnant black hole just after merger. The left panel 
of Fig. shows the horizon tendicity and the tendex 
lines emerging from the remnant black hole at some par- 
ticular time. The tendex structure is rotating counter- 
clockwise around the hole's horizon. The rotating tendex 
lines acquire accompanying vortex lines as they extend 
into the wave zone in a pattern like that shown in the 
right panel of Fig. and there they become gravita- 
tional waves. During merger, the horizon vortexes (right 
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FIG. 12: (Color online.) Remnant horizon, shown in the zy 
plane, just after merger, for a superkick simulation tuned for 
maximum remnant recoil in the +z direction; from a simu- 
lation presented in [20]. The black hole and vortex/tendex 
structures rotate counterclockwise. Left: Colors show hori- 
zon tendicity; tendex lines are shown emerging from the 
strongest horizon tendexes. Right: Colors show horizon vor- 
ticity; vortex lines are shown emerging from the strongest 
horizon vortexes. Figure adapted from [20]. See and 
for movies of horizon vortexes and tendexes for this simula- 
tion. 


panel of Fig. lock onto the same rotational angular 
velocity as the horizon tendexes, and generate gravita- 
tional waves in the same manner, with a pattern that 
looks like the left panel of Fig. 

In the wave zone, the gravitational waves produced 
by the rotating near-zone tendexes and those produced 
by the rotating near-zone vortexes superpose coherently, 
and the resulting radiation pattern depends on the angle 
between the horizon tendex labeled ^E" and the vortex 
labeled *B" in Fig. For the case shown, this angle is 
45 degrees, with Ex B in the —z direction (into the page). 
'This is the same as the structure of a gravitational wave 
propagating in the —z direction (Fig. [). As a result, 
the gravitational waves produced by the vortexes and 
those produced by the tendexes superpose constructively 
in the —z direction and destructively in the +z direction, 
resulting in a maximum gravitational-wave momentum 
flow in the —z direction and a maximum remnant recoil 
in the +z direction. 


4. Generic black-hole collisions 


The geometrodynamic behaviors of more general black 
hole collisions are now being explored numerically. For 
example, Fig. [13] shows trajectories of two black holes in 
a fully generic orbit. Vortexes from the larger, rapidly- 
spinning black hole pull the orbit of the smaller black 
hole into a complicated precession pattern. The spin di- 
rections of both black holes also precess as the holes or- 
bit each other. Eventually a common apparent horizon? 
forms around the two individual apparent horizons, and 
the two holes merge into one. The Ricci scalar (approxi- 
mately equal to -2 times the horizon tendicity) is shown 
on the two individual horizons and on the common hori- 
zon, at the moment the common horizon forms. 


5. Tidal disruption of a neutron star by a spinning black 
hole 


Our final example illustrates the interaction of ge- 
ometrodynamics and matter. Figure shows a simu- 
lation of a neutron star orbiting a black hole, a binary 
system important for gravitational-wave detectors and 
possibly for high-energy astrophysical phenomena such 
as gamma-ray bursts. Here the black hole has 3 times 
the mass of the neutron star, and a dimensionless spin 
S/M? = 0.5 in a direction inclined approximately 45 de- 
grees to the orbital angular momentum. When the orbit 


2 An apparent horizon is a surface of zero outgoing null expansion, 
and lies inside or on the event horizon. Apparent horizons are 
local quantities that are much easier to find in numerical simula- 
tions than are event horizons, because the location of the event 
horizon depends on the full future evolution of the spacetime. 


has shrunk sufficiently, because of energy lost to gravita- 
tional radiation, the black hole’s tidal tendexes rip apart 
the neutron star, and its frame-drag vortexes then pull 
the stellar debris out of its original orbit and into the 
black hole’s equatorial plane. If the black hole has a 
small enough spin or a large enough mass, the neutron 
star is not disrupted, but instead is swallowed whole by 
the black hole [32]. 


VI. GRAVITATIONAL-WAVE OBSERVATIONS 


Geometrodynamics generically produces gravitational 
waves. We are entering an era in which these waves, gen- 
erated by sources in the distant universe, will be detected 
on Earth. 

A first generation of interferometric gravitational-wave 
detectors has operated at sensitivities where it would re- 
quire luck to see waves. We were not lucky [35] [36]. 

Second generation detectors, much more advanced and 
complex in their design, are under construction. The first 
two of these (the advanced LIGO detectors in the U.S.) 
will begin operating in 2015 and should reach their de- 
sign sensitivity by 2019, and perhaps sooner [37]. They 
will be joined a bit later by the advanced VIRGO de- 


FIG. 13: (Color online.) The two thin curves are the trajecto- 
ries of the centers of two black holes in a generic orbit; from a 
simulation presented in [30]. The mass ratio of the two holes 
is 6:1, and the dimensionless spins of the larger and smaller 
holes are S/M? — 0.91 and 0.3, respectively (compared to 
the maximum possible spin S/M? — 1). The initial black 
hole positions are shown in black, and the initial spins are 
shown as arrows. The spins are initially oriented in generic 
directions, so that the orbital plane precesses. Shown also 
are the apparent horizons of both holes, and the common ap- 
parent horizon that encloses them, at the time the common 
apparent horizon first forms. The horizons are colored by the 
scalar Ricci curvature, which is approximately —2 times the 
horizon tendicity. See [31] for a movie of this simulation. 


(c) (d) 


FIG. 14: (Color online.) Four successive snapshots of a colli- 
sion between a black hole (black) and a neutron star (blue or 
gray), viewed edge-on in the initial orbital plane; from a sim- 
ulation reported in [33]. (a) The initial black-hole spin S is 
inclined with respect to the initial orbital angular momentum 
L. (b) The black hole's tendexes begin to rip apart the neu- 
tron star. (c) Some of the matter falls down the black hole, 
but some remains outside the horizon. (d) The black hole's 
vortexes pull the remaining matter into the black hole's equa- 
torial plane, forming a disk and a tidal tail. Figure adapted 
from [33]. See for a movie of this simulation. 


tector in Europe, KAGRA in Japan, and an advanced 
LIGO detector in India [38]. These advanced detectors 
will cover a 1000 times larger volume of the universe than 
the initial detectors did, with estimated event rates for 
mergers of black-hole and neutron-star binaries ranging 
from a few per year to a few per week [39]41]. Plans are 
being developed for further improvements, which should 
increase the event rate by another factor ten or more, 
without major changes in detector design. 

The most interesting gravitational-wave sources for 
these detectors, we think, are the dynamically evolving 
vortexes and tendexes attached to merging black holes, 
and to a black hole tearing apart a neutron star — the 
geometrodynamic phenomena discussed above. 

The numerical relativity community is building a cat- 
alog of binary simulations and associated gravitational 
waveforms to underpin the advanced detectors’ searches 
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for these waves. Simulations of binary black holes with 
several hundred different sets of parameters (mass ratios 
and initial vectorial spins) have been completed [30] [42]- 
[46]; and many more are underway or planned. A sample 
of computed waveforms from our SXS collaboration is 
shown in Fig. 

Once waves are being detected, comparison of observed 
waveforms with those from simulations will be crucial for 
understanding the waves’ sources. By such comparisons, 
can we deduce the sources’ geometrodynamics and test 
general relativity's predictions. 


VII. CONCLUSIONS 


Physicists have barely scratched the surface of ge- 
ometrodynamics. As numerical simulations continue to 
improve and are used to explore more complicated and 
generic situations, we expect to learn more about the 
geometrodynamics of critical behavior, singularities, dy- 
namical black holes, and other phenomena. We look for- 
ward to observations of gravitational waves from strongly 
gravitating astrophysical sources, which will enable us 
to test the geometrodynamical predictions of Einstein's 
equations for the first time. 
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FIG. 15: (Color online.) A recent catalog of black-hole binary simulations from [30]. Shown are 174 waveforms, each with two 
polarizations (shown in two colors or shades of gray) in a sky direction parallel to the initial orbital plane. The unit on the 
time axis corresponds to 0.1 s for a binary with a total mass of 20 solar masses. 
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